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<N ■ Abstract 

D ' In association with the Goldreich-Juhan's magnetic braking mechanism to explain the luminosity 
^ . 

. of radio pulsars as a result of their spin-down, it is of some interest to explore how much of magnetic 
flux can actually penetrate the surface of the rotating neutron stars at least in idealized situations. 

m 

^> ' In order to address this issue, one needs to figure out the amount of charge on and the structure 

OO . of magnetic field around a rotating neutron star. In the present work, based on the solution- 
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generating method given by Wald, the magnetic fields around both the uncharged and (slightly) 
charged neutron star have been obtained. Particularly for the charged neutron star, it has been 



• demonstrated following again the argument by Wald that the neutron star will gradually accrete the 

Q . charge until it reaches the equilibrium value Q = 2BqJ . Then next, the magnetic flux through one 
H , 

^ \ half of the surface of the rotating neutron star has been computed as well. With the nonvanishing 



accretion charge having value in the range < Q < 2BqJ, the total magnetic flux through the 



^ . neutron star has been shown to be greater than that without the accretion charge. 
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I. Introduction 



Radio pulsars are perhaps tlie oldest-known and the least energetic among all of the pul- 
sar categories. The thoretical study of these radio pulsars can be traced back to the 1969 
work of Goldreich and Julian ^ . In their pioneering work, Goldreich and Julian argued that 
pulsars, which are thought to be the rotating magnetized neutron stars, must have a mag- 
netosphere with charge-separated plasma. They then demonstrated that an electric force 
which is much stronger than the gravitational force will be set up along the magnetic field 
and as a result, the surface charge layer cannot be in dynamical equilibrium. Then there 
appear steady current flows along the magnetic field lines which are taken to be uniformly 
rotating since they are firmly rooted in the crystalline crust of the pulsar surface. Although 
it was rather implicit in their original work, this model for the pulsar electrodynamics in 
terms of the force-free magnetosphere does indeed suggest that the luminosity of radio pul- 
sars is due to the loss of their rotational energy (namely the spin-down) and its subsequent 
conversion into charged particle emission which eventually generates the radiation in the 
far zone. To elaborate on this last point, one should address the issue of what exactly the 
working mechanism for the release of the extracted rotational energy to infinitely far region 
would be. Indeed, if the field strength is large enough, say, of order B ~ lO^^(G') like the 
one observed near the radio pulsars, the vacuum surrounding the pulsar becomes unstable 
because any stray charged particles will be electrostatically accelerated and will radiate and 
this radiation will in turn produce further charged particles in the form of electron-positron 
pairs. When charges are produced so freely, the electromagnetic field around the neutron 
star will become approximately force-free, namely a force-free magnetosphere will be es- 
tablished. Therefore in reality, if the surrounding magnetic field is strong enough, charged 
particles are available due to the cascade production of electron-positron pairs and hence 
the extracted rotational energy of the neutron star may be released in the form of charged 
particle emission. Therefore in association with this Goldreich- Julian's mechanism for the 
spin-down of rotating neutron stars to implicitly explain the pulsar luminosity, it is of some 
significance to explore how much of the magnetic field lines actually penetrate the surface 
of the neutron stars at least in idealized situations. 

Upon appreciating the idea of Goldreich and Julian for pulsar magnetosphere, the next 
natural question would be ; what happens if we switch the central compact object from 
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rotating neutron stars to rotating (Kerr) black holes ? Will there be force-free black hole 
magnetosphere instead ? Indeed, such question has been asked and answered by Blandford 
and Znajek ^ long ago which since then has been known as the Blandford-Znajek mecha- 
nism. The way how some of the hole's rotational energy is carried off by the magnetic field 
lines threading the hole in the conventional Blandford-Znajek process can be summarized 
as follows. There the angular momentum and the rotational energy transport is achieved in 
terms of the conservation of "electromagnetic" angular momentum and energy flux flowing 
from the hole to its nearby field particularly when the force-free condition is satisfied in the 
surrounding magnetosphere. It is also amusing to note that the Blandford-Znajek mecha- 
nism can be translated into a simple-minded circuit analysis. In this alternative picture, 
the current, flowing from near the pole toward the equator of the rotating black hole geom- 
etry and the poloidal magnetic field together generate a "magnetic braking" torque which 
is antiparallel to the hole's spin direction. As a result, the rotating hole, which is a part of 
the circuit, spins down and its angular momentum and the rotational energy are extracted. 
For more details, we refer the reader to the original work of Blandford and Znajek and 
some comprehensive review articles 0] on this topic including the "membrane paradigm" 
approach 4]. Therefore, both the Goldreich- Julian and the Blandford-Znajek mechanisms 
can be dubbed "magnetic braking" mechanisms to explain the pulsar luminosity and the 
luminosity of active galactic nuclei (AGNs) respectively as a result of their spin-downs. As 
such, in order for both the Goldreich- Julian and the Blandford-Znajek processes to work, 
the local magnetic field lines should be firmly rooted on the surface of the central compact 
objects to carry away some of their rotational energies by means of this "contact interac- 
tion". 

In our earlier work P], it has been explored how much of the magnetic flux can actually 
penetrate the horizon of a rotating black hole in association with the Blandford-Znajek 
mechanism for rotational energy extraction from a Kerr hole. And in doing so, there we 
made use of the algorithm suggested by Wald to figure out the amount of charge on and 
the structure of magnetic field around a rotating black hole. The results concerning these 
issues presented there for the case of a rotating black hole can therefore be compared with 
those that will be presented here in this work for the case of a rotating neutron star since 
we shall in this work employ the algorithm of Wald's to determine the amount of accretion 
charge on and the structure of magnetic field around the rotating neutron star as well. Thus 
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it should be interesting to read the present work in parallel with our earlier work j5| on the 
Kerr black hole electrodynamics. 



II. The Hartle-Thorne metric for the region exterior to the slowly-rotating 

neutron stars 

We first begin by providing the rationale for treating the vicinity of a rotating neutron 
star relativistically as a non-trivial curved spacetime. Since the pioneering proposals by 
Gold y and by Pacini |2j, it has by now been widely accepted that indeed pulsars might 
be rotating magnetized neutron stars. Nevertheless, since then nearly all the studies on the 
pulsar electrodynamics have been performed by simply treating the region surrounding the 
rotating neutron stars as being flat. This simplification may be sufficient just to gain some 
insight into rough understanding of the origin of the pulsars' radiation. However, nearly 
over thirty years have passed of the study of pulsar electrodynamics and it seems to be 
that now we should look into the situation in a more careful and rigorous manner. In the 
present work, we would like to take one step forward in this direction. We are particularly 
aimed at exploring how much of the magnetic flux can actually penetrate the surface of a 
rotating neutron star and to this end we shall first determine the amount of charge on and 
the structure of magnetic field around a slowly-rotating neutron star by treating its exterior 
region literally as a curved spacetime and employing the algorithm of Wald's as mentioned 
in the introduction. Indeed, thus far only a handful of earlier studies concerning the elec- 
trodynamics of the rotating neutron stars appeared in the literature and they include the 
works by Muslimov and Tsygan, Muslimov and Harding and by Prasanna and Gupta 
j^. Although their motivation, approach and tools employed are completely different from 
those of our present work, RezzoUa et al. also considered analytic solutions of Maxwell 
equations in the internal and external regions of a slowly-rotating magnetized neutron star. 
Indeed, we now have a considerable amount of observed data for various species from 

f] n 

radio pulsars [l| to (anomalous) X-ray pulsars Even if we take the oldest-known 

radio pulsars for example, it is not hard to realize that these objects are compact 
enough which really should be treated in general relativistic manner. To be more con- 
crete, note that the values of the parameters characterizing typical radio pulsars are ; 
ro(radius) ~ lO^(cm), M(mass) ~ IAMq ~ 2x10^^(5'), r(pulsation period) ~ 10"^ — l(sec), 
i?(magnetic field strength) ~ 10^^ (G). Thus the Schwarzschild radius (gravitational radius) 



of a typical radio pulsar is estimated to be rsch = 2GMq/c^ ~ 3 x lO^(cm) (where G and 
c denote the Newton's constant and the speed of light respectively) and hence one ends 
up with the ratio ro/rsch ~ 10^(cm)/3 x lO^(cm) ~ 3. This simple argument indicates 
that indeed even the radio pulsar (which is perhaps the least energetic among all of its 
species) is a highly self-gravitating compact object that needs to be treated relativistically. 
Of course, the "relativistic" treatment here means that the region surrounding the pulsar, 
namely the pulsar magnetosphere has to be described by a curved spacetime rather than 
simply a flat one. Then the question now boils down to ; what would be the relevant 
metric to describe the vicinity of a rotating neutron star ? Although it does not seem to 
be well-known, fortunately we have such a metric of the region exterior to slowly-rotating 
relativistic stars such as neutron stars, white dwarfs and supermassive stars and it is the 
one constructed long ago by Hartle and Thorne Thus the Hartle-Thorne metric is a 
stationary axisymmetric solution to the vacuum Einstein equation and in the present work, 
we shall take this Hartle-Thorne metric as a relavant one to represent the spacetime exterior 
to slowly-rotating neutron stars. One might then wonder how the Hartle-Thorne metric 
for slowly-rotating neutron stars can be used to describe the millisecond pulsars which are 
thought to be rapidly-rotating magnetized neutron stars. Here, the "slowly-rotating" really 
means that the relativistic stars like neutron stars are relatively slowly rotating compared 
to the equal mass Kerr black hole which can rotate arbitrarily rapidly up to the maximal 
rotation J = M^. Thus as we shall discuss later on in the concluding remarks, the angular 
speed of the rotating neutron star represented by the Hartle-Thorne metric is indeed large 
enough to describe the millisecond pulsars. Now the Hartle-Thorne metric is given by 
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with M, J and Q being the mass, the angular momentum and the mass quadrupole moment 
of the (slowly) rotating neutron star respectively, P2{cos6) = [Scos^O — l)/2 being the 
Legendre poljTiomial, and being the associated Legendre polynomial, namely. 
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As is well-known, the only known exact metric solution exterior to a rotating object is the 
Kerr metric . Thus it would be worth clarifying the relation of the Hartle-Thorne metric 
br slowly-rotating relativistic stars given above to the Kerr metric. As Hartle and Thorne 
1^ pointed out, take the Kerr metric given in Boyer-Lindquist coordinate and expand it 
to second order in angular velocity (namely, the angular shift jS'^) followed by a coordinate 
transformation in the (r, 6')-sector, 
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where a = J/M. Then one can realize that the resulting expanded Kerr metric coincides 
with the particular case Q = J'^/M (with Q being the mass quadrupole moment of the ro- 
tating object) of the Hartle-Thorne metric. Therefore, in general this Hartle-Thorne metric 



is not a slow-rotation limit of Kerr metric. Rather, the slow-rotation limit of Kerr metric 
is a special case of this more general Hartle-Thorne metric. As a result, the Hartle-Thorne 
metric with an arbitrary value of the mass quadrupole moment Q can generally describe a 
(slowly-rotating) neutron star of any shape (as long as it retains the axisymmetry) . 
Now we turn to the choice of an orthonormal tetrad frame. And we shall particularly 
choose the Zero-Angular-Momentum-Observer (ZAMO) frame which is a sort of fidu- 
cial observer (FIDO) frame. Generally speaking, in order to represent a given background 
geometry, one needs to first choose a coordinate system in which the metric is to be given 
and next, in order to obtain physical components of a tensor (such as the electric and mag- 
netic field values), one has to select a tetrad frame (in a given coordinate system) to which 
the tensor components are to be projected. As is well-known, the orthonormal tetrad is 
a set of four mutually orthogonal unit vectors at each point in a given spacetime which 
give the directions of the four axes of locally-Minkowskian coordinate system. Such an or- 
thonormal tetrad associated with the Hartle-Thorne metric given above may be chosen as 
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The local, stationary observer at rest in this orthonormal tetrad frame has the worldline 
given by {dr = 0, dd = 0, {d<^ + fi^dC) = 0} which is orthogonal to spacelike hypersurfaces 
and has orbital angular velocity given by 



This is the long-known Lense-Thirring precession [15| angular velocity arising due to the 
"dragging of inertial frame" effect of a stationary axisymmetric spacetime. Indeed, it is 
straightforward to demonstrate that this orthonormal tetrad observer can be identified with 
a ZAMO carrying zero intrinsic angular momentum with it. To this end, recall that when 
a spacetime metric possesses a rotational (azimuthal) isometry, there exists associated ro- 
tational Killing field ip'^ = {d/dipY = ^ip such that the inner product of it with the tangent 
(velocity) vector = dx^/dr (with r denoting the particle's proper time) of the geodesic 
of a test particle is constant along the geodesic, i.e., 

L = gapi^'^uf' = g^t^'f'u' + g^^i^^u'*' 

dt d(j) , . 

= 9,tj^ + mj^- (10) 

Now, particularly when the local, stationary observer, which here is taken to be a test 
particle, carries zero angular momentum, L = 0, its angular velocity becomes 

d(j) {d(j)/ dr) gt^ 
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and this confirms the identification of the local observer at rest in this orthonormal tetrad 
frame given above with a ZAMO. 

III. Solution Generating Method by Wald 

1. Wald field 

From the general properties of Killing fields ; a Killing vector in a vacuum spacetime 
generates a solution of Maxwell's equations in the background of that vacuum spacetime, 
long ago, Wald [17] constructed a stationary axisymmetric solution of Maxwell's equations 
in Kerr black hole spacetime. To be a little more concrete, Wald's construction is based on 
the following two statements : 

(A) The axial Killing vector ip'^ = [d/dcj))^ generates a stationary axisymmetric test electro- 
magnetic field which asymptotically approaches a uniform magnetic field, has no magnetic 
monopole moment and has charge = 4J, 

= dip 

where "ci" denotes the exterior derivative and J is the angular momentum of a Kerr black 
hole. 

(B) The time translational Killing vector = {d/dt)^ generates a stationary axisymmetric 
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test electromagnetic field which vanishes asymptotically, has no magnetic monopole moment 
and has charge — —2M, 
F^ = dC 

where M is the mass of the Kerr hole. 

Apparently, this solution generating method of Wald's can be applied in a straightforward 
manner to the construction of a stationary axisymmetric solution of Maxwell's equations in 
a spacetime surrounding a slowly-rotating neutron star represented by the Hartle-Thorne 
metric given earlier. Thus we now look for the solution of the electromagnetic test field F 
which occurs when a stationary axisymmetric slowly-rotating neutron star of mass M and 
angular momentum J (represented by the Hartlc-Thorne metric) is placed in an originally 
uniform magnetic field of strength Bq aligned along the symmetry axis of the neutron star. 
And for now we consider the case when the electric charge is absent. Obviously, then the 
solution can be readily written down by referring to the statements (A) and (B) above as 
follows 
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and for the Hartle-Thorne metric given in the previous subsection, the solution above can 
be written in a concrete form as 
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Thus this can be thought of as the Wald-type field for the case of a rotating neutron star. Of 
course, a caution needs to be exercised and it is the point that this Maxwell field is correct 
only for a slowly-rotating neutron star but not for a neutron star with arbitrary rotation. 
2. Wald chcirge 

We now turn to the issue of charge accretion onto the (slowly-rotating) neutron star im- 
mersed in a magnetic field surrounded by an ionized interstellar medium ("plasma"). We 



10 



shall essentially follow the argument again given by Wald ^| and to do so, we first need 
to know the physical components of electric and magnetic fields. This can be achieved by 
projecting the Maxwell field tensor given above in eq.(14) onto the ZAMO tetrad frame. 
Thus using the dual to the ZAMO tetrad, = (cq = e(t),ei = 62 = 6(51), 63 = e((^)), 
given earlier in eq.(8), we can now read off the ZAMO tetrad components of F^^ as 



'S 



"20 — -^/xi^eaCo 
1 



rAV2 
F30 = F^^e^e'^ = 0, 

Fu = F^^ete'^ = 0, 

i^i3 = F.,ete^ 



2J_ 

M 



{de9<t>t 



(23) 



F23 = F^.e^e^ 



where, {dr9t^), {dr9tt), {de9t<t>), {de9tt), {dr9<t><i>) ^"^^ {de9^^) are as given above in eqs. (15)- 
(20). 

Here, consider particularly the radial component of the electric field (as observed by a local 
observer in this ZAMO tetrad frame) which, along the symmetry axis {6 = 0, it) of the 
slowly-rotating neutron star, becomes 
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Note that from this expression for the radial component of the electric field, it is not possible 
to determine the sign of charge accreted in terms of the relative direction between the angular 
momentum of the neutron star and the external naagnetic field. And this point is in contrast 
to what happens in the case of Kerr black hole ^, Q| studied previously. To be a little more 
concrete, there it has been realized that the radial component of the electric field for the case 
of Kerr black hole is radially inward/ outward if the hole's axis of rotation and the external 
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magnetic field are parallel/ antiparallel. Put differently, this implies that if the spin of the 
hole and the magnetic field are parallel, then positively charged particles on the symmetry 
axis of the hole will be pulled into the hole while if the spin of the hole and the magnetic 
field are antiparallel, negatively charged particles on the symmetry axis of the hole will be 
pulled into the hole. In this manner a rotating black hole will "selectively" accrete charged 
particles until it builds up "equilibrium" net charge. For the present case of rotating neutron 
star, however, an argument of this sort cannot be presented due essentially to the added 
complexity of the Hartle-Thorne metric compared to the Kerr metric. Nevertheless, one can 
still ask and answer the next natural question ; how then can the equilibrium net charge 
be determined ? And to answer this question, we resort to the "injection energy" argument 
originally proposed by Carter for the Kerr black hole case. Again, it can be applied 
equally well to the present case of (slowly-rotating) neutron star. Recall first that the energy 
of a charged particle in a stationary spacetime with the time translational isometry generated 
by the Killing field = (d/dt)'^ in the presence of a stationary electromagnetic field is given 
by 

with p^ = fhu^ — eA'^ being the 4-momentum of the charged particle with mass and charge 
m and e respectively. Now if we lower the charged particle down the symmetry axis into 
the rotating neutron star, the change in electrostatic energy of the particle will be 

^ final ^initial C^q,^ \horizon C^q,^ |oo 

Now, if < — > it will be energetically favorable for the neutron star to accrete particles 
with this charge whereas if > ^ it will accrete particles with opposite charge. In 
either case, the rotating neutron star will selectively accrete charges until A'^ is changed 
sufficiently that the electrostatic "injection energy" 6e is reduced to zero. We are then 
ready to determine, by this injection energy argument due to Carter, what the equilibrium 
net charge accreted onto the neutron star would be. In the discussion of Wald-type field 
given above, we only restricted ourselves to the case of solutions to Maxwell equation in the 
background of uncharged stationary axisymmetric neutron star spacetime and it was given by 
eq.(12). Now we need the solution when the stationary axisymmetric neutron star is slightly 
charged via charge accretion process described above. Then according to the statement (B) 
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in the discussion of Wald field given earlier, there can be at most one more perturbation of a 
stationary axisymmetric vacuum (neutron star) spacetime which corresponds to adding an 
electric charge Q to the neutron star and it is nothing but to linearly superpose the solution 
{-Q/2M)F(. = {-Q/2M)d^ to the solution given in eq.(12) to get 
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Then the electrostatic injection energy can be computed as 
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and ipfj,^'^ = (on the symmetry axis 6 = 0, tt) and tq denotes the radius of the neutron 
star along the symmetry axis. Thus one may conclude that a rotating neutron star in a 
uniform magnetic field will accrete charge until the gauge potential evolves to a value at 
which 6e = yielding the equilibrium net charge as Q = 2BqJ . Thus this amount of charge 
may be called Wald-type charge and it is particularly interesting to note that it turns out 
to be the same as the "Wald charge" (the equilibrium net charge) for the case of Kerr black 



hole 



17|. And this result appears to imply that the amount of equilibrium net charge 



Q = 2BqJ might be a generic value once a general relativistic rotating object gets charged 
via accretion. There, however, does exist a distinction between the case of Kerr black hole 
and that of a rotating neutron star. That is, from eq.(28), the rotating neutron star will 
eventually accrete the equilibrium net charge oi Q = 2BqJ only if [1 — Ai?|r=ro,e=o,7r] 7^ 0. 
Certainly, this feature was absent in the case of Kerr black hole and can be attributed to the 
different structure of Hartle-Thorne metric describing the region surrounding the rotating 
neutron star. Therefore, particularly if the values of mass (M), angular momentum (J), 
mass quadrupole moment (Q) and radius along the symmetry axis (rg) are such that they 
satisfy 
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M3 
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then it is "energetically unlikely" that a slowly-rotating neutron star would accrete any 
charge and we find it a peculiar feature which has no parallel in the Kerr black hole case. It, 
however, should be taken with some caution (i.e., should not be taken as being conclusive) 
since the injection energy argument due to Carter that leads to this conclusion is, rigorously 
speaking, indeed restricted to the situation along the symmetry axis 6 = 0,tt and at this 
point it is not so obvious that the same would be true off the symmetry well. 
Next, recall that we announced from the beginning that we shall consider the case when the 
charge accreted on the neutron star is small enough not to distort the background Hartle- 
Thorne geometry. Now we provide the rationale that this is indeed what can actually 
happen. To do so, we first assume that just like in the case of Kerr black hole (a particular 
slow-rotation limit of which is the Hartle-Thorne geometry), in the present Hartle-Thorne 
spacetime, J < M^. Then using the fact that the typical value of the charge accreted on the 
rotating neutron star is Q = 2BqJ as just described, its charge-to-mass ratio has an upper 
bound 

where in the last equality we have converted Bn and M from geometrized units to solar-mass 

n 

units and gauss |5]. Thus for the typical case of radio pulsar ; with mass M ~ Mq in the 
surrounding magnetic field of strength Bq ~ lO^^(G'), Q /M ~ 10^^ << 1. Thus as we can 
see in this radio pulsar case, the charge-to-mass ratio of the associated rotating neutron star 
is small enough not to disturb the geometry itself. Thus we can safely employ the solution- 
generating method suggested by Wald to construct the solution to the Maxwell equations 
when some amounts of charges are around to which we now turn. 

IV. Stationary axisymmetric Maxwell field around a "slightly charged" and 

slowly-rotating neutron star 

As discussed in the previous subsection, the stationary axisymmetric solution to the 
Maxwell equation in the background of a rotating neutron star with charge Q accreted in a 
originally uniform magnetic field can be constructed as 



14 



Again for the Hartle-Thorne metric for a slowly-rotating neutron star given in the previous 
subsection, the solution above can be written as 



and here, {drgt^p), {drOtt), {degt^), {degtt), {drg<i><i>) and {deg^,^,) are as given earlier in the 
previous section. 

Obviously, in order to have some insight into the nature of this solution to the Maxwell 
equations, one may wish to obtain physical components of electric field and magnetic in- 
duction. And this can only be achieved by projecting the Maxwell field tensor above onto 
an appropriate tetrad frame like ZAMO tetrad which, for the Hartle-Thorne metric being 
considered, can be well-defined as mentioned earlier. The ZAMO tetrad is well-known and 
widely employed in various anaysis in the literature. ZAMO is a fiducial observer follow- 
ing a timelike geodesic orthogonal to spacelike hypersurfaces. And this implies that its 
4- velocity is just the timelike ZAMO tetrad u*^ — Cq given earlier in eq.(8). Thus upon 
projecting the Maxwell field tensor components on the ZAMO tetrad frame, the physical 
components of electric and magnetic fields can be read off as Ei = Eiq = Efj_i^{efeQ) and 
Bi = = eijkFf"'{e'^el)/2, respectively. Namely, from the dual to the ZAMO tetrad, 

eA = (eo = e(t), ei = e(^r), 62 = e(e), 63 = 6(0)), given earlier in eq.(8). 




(31) 



10 




(|)"'|b»-"!i!.{o,..i»|(i- 




) 



) 



(drgtt) 
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i^30 = F^^e^e'^ = 0, (32) 

Fu - i^M-e^e^ = 0, 

-^13 = F^iyCie'^ 

Ax^/^ 1 1 f 2J f Q \ ] 



F23 = F^.e^e^ 



V. The magnetic flux through the neutron stcir 

With the asymptotically uniform stationary axisymmetric magnetic field which is aligned 
with the spin axis of a "slightly charged" slowly-rotating neutron star given above, we now 
would like to compute the flux of the magnetic field across one half of the surface of the 
neutron star which is assumed to be of exact spherical geometry with radius tq. The physical 
motivation behind this study is to have some insight into the question of how much of the 
magnetic flux can actually penetrate the surface of a rotating neutron star - at least in 
idealized situations. We now begin by considering two vectors lying on the surface of the 
neutron star which are given by 

dx'^ = (0, 0, d9, 0), dx"^ = (0, 0, 0, dcp). (33) 

Then in terms of the 2nd rank tensor constructed from these two vectors [10], 

d(7°^ = ^ {dx^x^ - dx(dx^) (34) 

one now can define the invariant surface element of the neutron star as 

ds = (2d(7«^a"^)i/2 (35) 
= {geegHf^ Iro Md(t>- 

Next, since the tensor da"'^ is associated with the invariant surface element of any (not 
necessarily closed) 2-surface, the fiux of electric field and magnetic field across any 2-surfaces 
can be given respectively by 

= j F^pda^P = I ]^e^p^'F,sda^P, 

^ I F^f.da^f'. (36) 
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In particular, the flux of magnetic fleld across one half of the surface of the neutron star is 



al3 



2n 



tt/2 



(37) 



Then using F^^ evaluated on the surface of the neutron star 

Q ^ 



'-Bo 
2 " 



2J / 



2^0 J 



(38) 



Q 



4J 



M 

5 Q- J^/M 
8 M3 



+ 



2BoJ ) ro 
2M 



1 + 



J2 



Mrl 



1 + 



2M\ 



[r2(l-2M/ro)]i/2^2(^ 1) 1) 



sm y cos I 



+ l^Bo{2rt 



2_ 2J 



4J 



^ 5 Q - JVM I 



2BoJ; ro 
2M 



J2 



1 + ?^) 



8 M3 



\+[r2(l-2M/ro)]i/2^2(5 1) QK^ 1)^ 



(sin cos^ 6 — siv? 9 cos ( 



and the result of integration, ^^^"^ dO sin 6 cos 6 = 1/2, /g'^^^ (i^^(sin cos'^ 6* — sin'^^cosl 
we flnally arrive at 



$B = -kBo Y'o 



1 



Q 



Mrg 



Mro \ 2BoJ ^ 
2M\ 5Q-J^/M 



X 



(39) 



8 M3 



+ 



2M 



T7J0J(S - 1) - - 1) 



M 



[r2(l - 2M/ro)]i/2 

This is the magnetic flux through one half of the surface of the slowly-rotating neutron star 
with accretion electric charge Q. Some discussions concerning the nature of this magnetic 
flux are now in order : 



(i) In the absence of the accretion charge, i.e., Q = 0, the magnetic flux above appears 
to drop to zero (i.e., the magnetic fleld gets entirely expelled) when the neutron star is 
maximally-rotating, which, based on the expression for the magnetic flux across the surface 
of the neutron star given in eq.(39), amounts to J^„^ = Mtq/A. This curious feature 
indeed has its counterpart in the Kerr black hole case |l9| which has been known for some 
time. For the present case of neutron star, however, we are not fully qualifled to make 
any deflnite statement on this point since the Hartle-Thorne metric employed represents 
the region surrounding a slowly-rotating neutron star. Thus the expression for the mag- 
netic flux given above may be invalidated for the case of rapidly-rotating neutron star and 
that is beyond the scope of the present work as the associated metric is not available for now. 
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(ii) With the nonvanishing accretion charge having value in the range < Q < 2BoJ, it 
is interesting to note that 

^ 0) > $b(Q = 0). (40) 

Actually, this is one of the points of central importance we would like to make in the present 
work. The physical interpretation of this characteristic can be briefly stated as follows. 
Once the neutron star accretes charge, toroidal currents occur as a result of rotation of the 
star. This toroidal currents, in turn, generates new magnetic fields which would be additive 
to the existing ones increasing the total magnetic flux through the surface of the rotating 
neutron star. As has been studied in detail in Q recently, essentially the same phenomenon 
takes place in the case of slightly charged Kerr black hole. There, the interpretation was 
even clearer ; when the spin of the hole and the asymtotically uniform magnetic field are 
parallel (antiparallel), the hole selectively accretes positive (negative) charges according to 
the injection energy argument proposed by Carter and they, in turn, generate magnetic 
fields which are clearly additive to the existing ones. And we expect that essentially the 
same will be true for the present rotating neutron star case although we could not determine 
the sign of accretion charge in association with the relative direction between the spin of 
the neutron star and the direction of the magnetic field due to the highly complex structure 
of Hartle-Thorne metric as we mentioned earlier. 



(iii) For the Wald charge, Q = 2BqJ, the magnetic flux across the surface of the rotating 
neutron star appears to take on rather special value 

$B = BoTrrI X (41) 

As mentioned in the introduction, it may be interesting to compare this analysis of the 
nature of the magnetic flux through a rotating neutron star with that of the Kerr hole 
studied in P]. 

VI. Concluding remarks 

In the present work, based on the solution-generating method given by Wald, the magnetic 
fields around both the uncharged and (slightly) charged neutron star have been obtained. 
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Particularly for the charged neutron star, it has been demonstrated following again the ar- 
gument by Wald that the neutron star will gradually accrete the charge until it reaches the 
equilibrium value Q = 2BqJ . In association with the "magnetic braking" model of Gol- 
dreich and Julian l| for the spin-down of magnetized rotating neutron stars, the magnetic 
flux through one half of the surface of the rotating neutron star has been computed as well. 
With the nonvanishing accretion charge having value in the range Q < Q < 2BoJ, the total 
magnetic flux through the neutron star has been shown to be greater than that without 
the accretion charge. Indeed, the physical interpretation of this characteristic can be briefly 
stated as follows. When the spin of the neutron star and the direction of the asymtotically 
uniform magnetic field are parallel (antiparallel), the neutron star would selectively accrete 
positive (negative) charges following the injection energy argument proposed by Carter and 
they, in turn, generate magnetic fields additive to the existing ones. 

Next, one might be worried about the validity of the Hartle-Thorne metric for the region 
surrounding the slowly-rotating neutron stars employed in this work to describe the mag- 
netosphere of pulsars which seem rapidly-rotating having typically millisecond pulsation 
periods. Here the "slowly-rotating" means that the neutron star rotates relatively slowly 
compared to the equal mass Kerr black hole which can rotate arbitrarily rapidly up to the 
maximal rotation J = M^. Thus this does not necessarily mean that the Hartle-Thorne 
metric for slowly-rotating neutron stars cannot properly describe the millisecond pulsars. 
To see this, note that according to the Hartle-Thorne metric, the angular speed of a rotat- 
ing neutron star is given by the Lense-Thirring precession angular velocity in eq.(9) at the 
surface of the neutron star, which, restoring the fundamental constants to get back to the 
gaussian unit, is 

a; = ^(|), with J = aM^(^) (0 < a < 1). (42) 

Thus the Hartle-Thorne metric gives the angular speed of a rotating neutron star, hav- 
ing the data of a typical radio pulsar, M ~ 2 x 10^'^(5'), tq ~ 10^ (cm), as = 
2a{M'^ / rQ){G / c){G / c^) ~ 10'^(l/sec) which, in turn, yields the rotation period of r = 
27r/cj ~ 10-2(sec). And here we used, (G/c^) = 0.7425 x {cm / g) and (G/c) = 
2.226 X [crn^ / g ■ sec). Indeed, this is impressively comparable to the observed pul- 
sation periods of radio pulsars r ~ 10^^ — l(sec) we discussed earlier. As a result, we expect 
that the Hartle-Thorne metric is well-qualified to describe the geometries of millisecond pul- 
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sars. 

Lastly, we mention an inherent drawback of our treatment of the issues presented in this 
work. In the present work, we have mainly considered the case with symmetric geometry in 
which the stationary axisymmetric magnetic field is precisely aligned with the neutron star's 
axis of rotation. Although this simplification/idealization was inevitable in order to make 
the system mathematically tractable, this would render the magnetized rotating neutron 
star fail to describe pulsars in the real world as the pulsar involves, by definition, oblique 
spin-magnetic field configuration. (In the case involving a rotating black hole, however, such 
oblique configuration cannot be maintained since any possible initial misalignment between 
the spin and the direction of the magnetic field will be removed in an efficient manner either 
by a viscous torque when the accretion disc is involved [2^ or by a magnetic alignment 
torque j2l|.) Indeed, almost all of the theoretical studies of pulsar electrodynamics assume 
this idealized symmetric geometry for concrete, analytical calculations and concerning this 
limitation, our work makes no progress in this direction. Thus we are still away from the 
satisfactory theoretical description of the real physics that takes place in the pulsar mag- 
netosphere. Nevertheless it is our hope that the result of the present work, obtained in a 
more realistic situation when some amount of plasma (charges) are around in addition to the 
magnetic field, eventually lends supports to the operational nature of the Goldreich- Julian 
mechanism of the magnetic braking to explain the pulsar radiation due to its spin-down and 
puts it on a firmer ground at least on the theoretical side. 
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